Accuracy of the quantum capacitor as a single-electron source 
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A periodically driven quantum capacitor may function as an on-demand single electron source as it 
has recently been demonstrated experimentally. However, the accuracy at which single electrons are 
emitted is not yet understood. Here we consider a conceptually simple model of a quantum capacitor 
and find analytically the noise spectrum as well as the counting statistics of emitted electrons. We 
find that the failure rate of the capacitor can be arbitrarily small when operated under favorable 
conditions. Our theoretical predictions may be tested in future experiments. 

PACS numbers: 73.23.-b, 72.10.-d, 72.70.+m 
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Introduction. — Controllable single electron sources are 
at the forefront of current research on nano-scale elec- 
tronics. Systems that generate quantized electrical cur- 
rents, for example quantum capacitors [1:^3:] and quan- 
tum pumps fij , are of great interest due to their poten- 
tial applications in metrology [5] and quantum informa- 
tion processing [6] as well as in basic research on single- 
and few-electron physics in mesoscopic structures. The 
quantum capacitor constitutes one archetype of a single 
electron emitter. The capacitor consists of a nano-scale 
cavity that exchanges particles with a reservoir through 
a narrow constriction, a so-called quantum point contact. 
When the capacitor is subject to periodic voltage mod- 
ulations, single electron emission and absorption occur 
at giga-hertz frequencies as it has recently been demon- 
strated experimentally [1]. It has also been verified 
that the relaxation resistance of the capacitor is quan- 
tized in units of hjle^ independently of microscopic de- 
tails, in agreement with theoretical predictions [8|. The 
experiments indicate that the quantum capacitor may 
function as an on-demand electron source that ideally 
can be controlled down to the level of single electrons. 

Despite the experimental and theoretical advances, the 
accuracy at which the quantum capacitor emits electrons 
is still not well understood. This is an important question 
for potential applications and it may ultimately be the 
criterion that determines if the quantum capacitor be- 
comes an integrated part of future nano-scale electronics 
operating at giga-hertz frequencies. In this work we an- 
alyze the accuracy of the quantum capacitor as a single 
electron source. To this end, it is necessary not only to 
study the mean current of electrons emitted by the ca- 
pacitor but also the fluctuations of the current around 
the mean. We describe the capacitor using a simple 
model which has been shown to reproduce noise measure- 
ments in numerical simulations of the recent experiment 
reported in Ref. Q . We derive an analytic expression for 
the noise spectrum that fully accounts for the measure- 
ments. Furthermore, we characterize fluctuations in the 
current of emitted electrons by evaluating the counting 
statistics and find that the failure rate of the device under 
optimal operating conditions can be vanishingly small. 
































— ► 




FIG. 1: Quantum capacitor, a) An edge state is connected 
to the capacitor via a quantum point contact. The periodic 
voltage Uif) applied to the capacitor causes emission and ab- 
sorption of single electrons to and from the edge state, b) 
Periodic voltage ?7(i) and resulting current {I{t)) as func- 
tions of time. Numerical and analytical results are shown, 
c) Analytic result for the noise spectrum Vi{iJ — 2n /T) as 
function of the correlation time t (see text). Experimental 
results have been adapted from Ref. 0]. 



Quantum capacitor. — The system is shown schemat- 
ically in Fig. [T^. An edge state is connected to the ca- 
pacitor via a quantum point contact whose transmission 
probability p can be controlled with external gates. With 
the quantum point contact pinched off, the capacitor has 
a discrete energy spectrum with typical level spacing A, 
which is much larger than temperature, and Tq = /i/A is 
the time it takes an electron to travel a full round along 
the edge of the capacitor. As the quantum point contact 
is opened, electronic states of the capacitor with ener- 
gies below the Fermi level of the external reservoir are 
filled. We now consider the situation where a step-like 
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gate voltage U{t) periodically shifts the highest occupied 
level above and below the Fermi level of the reservoir. 
The period T is much longer than Tq and the amplitude 
2Uo is on the order of the level spacing A/e. This causes 
periodic emission of coherent single electron wave packets 
from the capacitor to the outgoing edge state, followed by 
refilling from the incoming edge state as it was recently 
demonstrated experimentally p|. In this work, we study 
the accuracy of single electron emission as the ratio of 
the period T and the escape, or correlation, time r (de- 
fined below) is varied or, equivalently, as function of the 
dimensionless parameter e = e"-^/^"^. 

Model. — The charge dynamics of the capacitor follows 
a simple model that recently was shown to reproduce 
measured data in numerical simulations of the experi- 
ment reported in Ref. . The absorption phase of dura- 
tion T/2 (denoted by ® in Fig.[TjD) is discretized in time 
steps of length Tq during each of which a single electron 
can enter the capacitor with probability p. The emission 
phase (denoted by (2)) is similarly discretized in time, and 
in each time step the probability of emitting an electron 
is p. With the amplitude of the periodic driving being 
on the order of the level spacing, higher-lying states can 
safely be neglected and only a single (additional) electron 
can occupy the capacitor. This semi-classical picture can 
be formulated as a master equation in discrete time for 
the probability of the capacitor to be occupied by an 
electron. Setting the electron charge e = 1 in the follow- 
ing, this probability is equal to the average (additional) 
charge of the capacitor (Q), where Q = 0, 1. The master 
equation determines the evolution of the average charge 
after one time step and reads 9] 

r p[i-{Qitk.i))]+{Qitk.i)) ® 

{Q{tk+i,i)) = { (1) 

[ {l-p){Q{tkA)) ® 

where we have used that 1 — (Q) is the probability for 
the capacitor to be empty and t = tk.i denotes time at 
the fc'th time step during the Vth period. The absorption 
(emission) phase ® ((2)) corresponds to k — 1,2...,K 
{K + 1,K + 2,..., 2K). Although, we do not derive 
the master equation here, some insights into its origin 
can be found by considering the current as it was cal- 
culated in Ref. [Ij using scattering matrices. The cur- 
rent was shown to consist of one step-like term with step 
length To and one oscillatory part with period Tq. The 
oscillatory part is due to quantum interference and van- 
ishes with increasing temperatures. At arbitrary temper- 
ature, only the first term remains after integration of the 
current over the time step Tq, which leads to the mas- 
ter equation above. As in the experiment [§] (with the 
measurement frequency equal to the driving frequency 
27r/a; — T = 60 Tq), we consider time scales that are 
much longer than Tq for which the master equation pro- 
vides a quasi-continuous description. 

Average charge and current. — It is straightforward to 



solve the master equation ([T]) for the average charge (Q) 
and thus the net current running out of the capacitor 
{I{t)) = -{Q{t)) ~ [{Qit))-{Q{t + To))]/To. The average 
charge can be cast in a form similar to that of an RC 
circuit reading 

r l-/3,e-(*'=-'-'^)/" (f) 

{QM = { (2) 

i a;e-(*^.<-['+^m/- (2) 

where we have defined ai — 1/(1 -I- e) + and /?/ — 
1/(1 + e) — Oe^^^^ with depending on the initial con- 
ditions at the time when the periodic driving is turned 
on. The correlation time r = To/ln[l/(l — p)] deter- 
mines the time scale over which the system loses mem- 
ory about the initial conditions encoded in 6 and {Q) 
becomes periodic in time. Figure [T}d illustrates the excel- 
lent agreement between the current {I{t)) obtained from 
the analytic expression ^ and our numerical simulations 
based on Eq. ([Ij, including the initial transient behav- 
ior. The mean charge emitted during the emission phase 
is tanh(r/4r) [see also Eqs. © and ©]. 

Noise spectrum. — In the experiment, the Fourier 
transform of the time-averaged correlation function 

{SI{t)6I{t + to))* with Slit) = I{t) - {I{t)) was mea- 
sured [9]. The noise spectrum is then Vi{lo) — 

dto {5I{t)5I{t + to)) *e'"*° or Viiuj) ~ uj^Vgiu) in 
terms of the corresponding charge correlation function 
Vq{lo). We evaluate the charge correlation function 
by noting that {Q{t)Q{t + to)) is the joint probabil- 
ity for the capacitor to be charged with one electron 
at time t and at time t + to. Using conditional prob- 
abilities jve write {Qit)Q(t + to)) = {Q{t)){Q{t + to)), 
where {Q{t + to)) is the probability that the capacitor 
is charged with one electron at time t + to given that 
it is charged at time t. For > 0, {Q{t + to)) can 
be found by propagating forward in time the condition 
{Q{t)) = 1 using the master equation ([T]). Similar rea- 
soning applies to the case to < 0. Integrating over t, 
the time-averaged charge correlation function becomes 

{dQ{t)dQ{t + to))* = ^e~l*''l/^tanh(^), and finally, 
we obtain the noise spectrum 

^^H-|tanh(|)^^. (3) 

Figure [TJ; shows our analytic expression for the noise 
spectrum together with experimental results adapted 
from Ref. [^. The analytic result captures the experi- 
ment over the full range of correlation times r and inter- 
polates between the two limiting cases discussed in Ref. 
^. In the shot noise regime r ^ T (e ~ 1), the prob- 
ability of emitting and reabsorbing an electron during a 
period is very small, and electron emission becomes rare. 
In this regime, we find Vi{uj) — > 1/2t in agreement with 
Ref. [9|]. In the phase noise regime r ^ T (e ~ 0), the 
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probability of emitting and absorbing an electron dur- 
ing each period is close to one, and the main source of 
finite- frequency fluctuations is the random times of emis- 
sion and absorption within a period. In this regime we 
find Viiuj) I ifjil'i as suggested in Ref. Q. The 
zero- frequency limit P/(0) — reflects that charge does 
not accumulate on the capacitor over time once (Q) has 
become periodic in time. 

Counting statistics. — While the noise spectrum ([3]) is 
related to the net current running out of the capacitor 
(including both absorption and emission of electrons) , it 
is relevant to characterize the emission process alone in 
order to quantify the accuracy of the capacitor as a sin- 
gle electron source. To this end, we consider the count- 
ing statistics of emitted electrons and hence introduce 
the probabilities Pj{n,l) that the capacitor is occupied 
by j = 0,1 (additional) electrons, while n electrons have 
been emitted after I periods. The probability distribu- 
tion for the number of emitted electrons is P{n, I) — 
Po(n, l)+Pi{n, I) and the corresponding cumulants ((n™)) 
of the distribution are defined through the cumulant gen- 
erating function (CGF) 5(x, I) = In P{n, /)e™x] as 
((n™)) = dJ^^^S{xJ)\x^o- The CGF can be written 
Six, I) = In [1 • P(x, I)], where 1 - [1, 1]^ and P(x, I) = 
[Piix,l),PoixJ)f with P,ix,l) = E„^:,KOe"^ J - 
0, 1. The evolution of the probability vector after one 
period of the driving is obtained from the master equa- 
tion (P) and reads P{x,l 
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, , and we have introduced the count- 

pe"^ 1 

ing field x in the off-diagonal element of L2 correspond- 
ing to electron emission [l^]. The CGF after I periods 
is then S{x,l) = In [l • A'(x)Pin] with Pjn being the 
X-independent initial condition as counting begins. For 
a large number of periods, A'(x) is dominated by the 
largest eigenvalue of A(x) 



A(x) =e + e'- 



1 



-IX 



(4) 



such that 5(x, ^ 'ln[A(x)] ^H]- The expression for 
the CGF is a powerful result that allows us to fully char- 
acterize fluctuations in the current of emitted electrons. 

Cumulants. — The cumulants of the current are defined 
as the constant ratio ((/'")) = after a large num- 

ber of periods. The first three cumulants are 

1 - i 



((/» 



1 



in) = 



2e 



+ 

2e(4£ 



tanh (T/4t) 



(5) 
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Higher-order cumulants can be approximated by noting 
that the CGF has square-root branch points at ix± = 



ln[(l — e)'^/Ae] ± in, close to which the CGF behaves as 
S{x, I) — 2Z-\/ix± — ix- Following Ref. [12] we find for 



large orders ((/™)) 



4B„ 



cos [(m- l/2)arg(ix-(-)] 



with Bfn,^ ~ /i(/i + 1) . . . (/i + m — 1). The asymptotic 
expression gives excellent agreement with exact results 
for the high-order cumulants as we have checked. 

Figure [2] illustrates the good agreement between nu- 
merical simulations of the first four cumulants and our 
analytical results. The first cumulant (the mean current) 
is equal to the average charge emitted during one period 
which varies from in the shot noise regime (e ~ 1) 
to 1 in the phase noise regime (e ~ 0). In the shot 
noise regime the first few cumulants can be written as 
((/")) ~ (1/2)™-! ((/)), indicating that the transport 
statistics is similar to that of a Poisson process with an 
effective charge of 1/2 - this is characteristic for a Pois- 
son process in which only every second event results in 
emission. In the phase noise regime, the mean current is 
close to 1 and the first few cumulants are close to zero 
since electrons arc emitted in an orderly manner due to 
the periodic driving. Low-frequency fiuctuations in the 
stream of emitted electrons arise only in the rare cases 
when the capacitor is not charged in the absorption phase 
or when it fails to emit in the emission phase. We pre- 
dict the occurrence of such failures in the phase noise 
regime by expanding the CGF to lowest order in e as 
S{x, I) — ^ix + 2e(e-'-*^ — 1)]. The first term corresponds 
to a deterministic process in which one electron is emit- 
ted in each cycle. The second term is the sum of two in- 
dependent Poisson processes describing "cycle missing" 
events occurring with rate e per period, either because 
the capacitor is not charged or because it fails to emit. Of 
course, if the capacitor is not charged it also fails to emit, 
however, such correlations do not enter to lowest order 
in e. In the cross-over from the shot noise to the phase 
noise regime, the third cumulant changes from positive 
to negative before it vanishes in the phase noise regime. 
The negative third cumulant signals a left-skewed dis- 
tribution caused by the mean current ((/)) being close 
to the upper limit of 1, which cuts off the distribution 
to the right. The fourth cumulant goes through a re- 
gion with negative values followed by a region with pos- 
itive values before vanishing in the phase noise regime. 
The negative (positive) fourth cumulant indicates a sub- 
gaussian (super-gaussian) distribution with light (heavy) 
tails. Higher-order cumulants become increasingly oscil- 
lating functions of e as recently predicted .12)]. 

Large deviation function. — The above statements can 
be further corroborated by calculating the full distribu- 
tion function P{I,l) J^^ ^e^{HHx)]-^xi} ^^le cur- 
rent / = n/l measured after many periods using a saddle- 
point approximation [lol [Tlj . This procedure yields 
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FIG. 2: Mean current and cumulants. Mean current 
((/)) together with the normalized cumulants of the current 
((/'")) /((/)), m = 2,3,4, as functions of the correlation time 



T or, equivalently, e 



-T/2t 



Symbols indicate results ob- 



tained from numerical simulations (obtained from 50.000 re- 
alizations), while dashed lines show our analytic expressions. 



FIG. 3: Large deviation function. The probability distribu- 
tion P(/, I = 50) of the current I for different values of the 
correlation time r or, equivalently, e = e"^^^'^. The full lines 
show the saddle-point approximation ((6]), while dashed lines 
correspond to a Gauss distribution with mean and variance 
(the second cumulant) obtained from Eq. The histograms 
were assembled from 50.000 numerical realizations. 



In Fig. 13] we show our analytic result together with nu- 
merical simulations. For comparison, we also show a 
Gauss distribution with mean and variance (second cu- 
mulant) obtained from Eq. ([5]). The analytic result de- 
scribes our numerics very well. Close to the phase noise 
regime (e ~ 0), the distribution function is clearly left- 
skewed in accordance with the third cumulant being neg- 
ative in Fig. [2j Additionally, the distribution is super- 
gaussian with heavy tails in comparison to the Gauss 
distribution. As the shot noise regime is approached, 
the distribution becomes sub-gaussian and light-tailed. 
These qualitative changes are reflected in the sign-change 
of the fourth cumulant seen in Fig. [2j In the phase noise 
regime, the distribution is 1 for 1 = 1 and zero oth- 
erwise. This can be understood by considering the ex- 
pected mean number of "cycle missing" events in our nu- 
merical simulations. After I = 50 periods of the driving, 
the expected mean number 2^6"^/^"^ is still vanishingly 
small for r/T = 0.01, and we do not observe a single "cy- 
cle missing" event during 50.000 numerical realizations as 
seen in Fig. [S] The large deviation function was recently 
measured in single electron transport through a Coulomb 
blockade quantum dot [l3|, and the results presented in 
Fig. |3]may serve as an experimental test of our model in 
similar measurements on a quantum capacitor. 

Conclusions. — We have analyzed the quantum capaci- 
tor and found analytically the noise spectrum which fully 
accounts for recent measurements. We have character- 
ized the accuracy of the capacitor as a single electron 
source through calculations of the counting statistics and 
found that the failure rate can be arbitrarily small under 
favorable operating conditions. Our results are impor- 
tant for possible applications of quantum capacitors in 
nano-scale electronics operating at giga-hertz frequencies 



and our predictions may be tested in future experiments. 
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